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Abstract 

The global in time existence of solutions of a system describing 
the interaction of gravitationally attracting particles with a general 
diffusion term and fixed energy is proved. The presented theory covers 
the case of the model with diffusion that obeys Fermi-Dirac statistics. 
Some of the results apply to the dissipative polytropic case as well. 
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1 Introduction 

We consider the following initial-boundary value problem 

rit = V ■ {D{Vp{n) +nV<^)) in x (0, oo) , (1.1) 

A(f = n in Qx (0, oo) , (1.2) 

(Vp(n) + Vv^) • z> = = on dQx{0,oo), (1.3) 

n{0)=no>0 in QcR'^, (1.4) 
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with the pressure in the self-similar form 



(1.5) 



for a given function P and some diffusion coefficient motivated by applica- 
tions to statistical mechanics and describing self-attracting clouds of particles 
modelling elliptical gallaxies, globular clusters, interstellar medium or cores 
of neutron stars among others (cf. [12] and references therein). These sorts of 
problems were considered among others in [S] including numerous pressure 
formulae coming from statistical mechanics: Maxwell-Boltzmann, Fermi- 
Dirac, Bose-Einstein and polytropic distributions. The common feature of 
all these examples is the self-similar profile of the pressure (11.51) . In this pa- 
per we focus our attention on the Fermi-Dirac model although we formulate 
the results in a more general setting. The pressure in this model assumes 
an intermediate form between the linear Maxwell-Boltzmann case at zero 
and a polytropic, power-like form at infinity. In [3] the authors proved the 
local and global existence for the specific choice of the diffusion parameter D 
corresponding to the Fermi-Dirac statistics in the isothermal case, i.e., with 
fixed, constant temperature 9. Moreover, the asymptotic behaviour with the 
possibility of the evolution towards steady states was adressed therein (cf. 
also [23j), whereas in [4j some results for nonisothermal case were established. 
In the aforementioned papers, some a priori estimates for the density n and 
the pressure P were also provided. For physical motivations one can see 
the series of papers of Chavanis and collaborators including [7]-[l5]. Related 
Keller-Segel model in mathematical biology was recently studied in [18], [26] 
and [27] and the blow-up for large data was proved. 
Note that, as a consequence of (11.31) . total mass 



is conserved during the evolution of the system. 

In the first part of the paper we will extend the results of [3] to allow 
more general pressure and nonconstant temperature for the dimension d G 
[2, 4] with small mass M. The local and global results for the parabolic 
perturbation is contained in Proposition 12. H and the result for the original 




(1.6) 
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parabolic-elliptic problem can be found in Theorem 13.11 As far as steady 
states are concerned nonexistence results hold for d > 2(1 + y^) (cf. [21j) we 
can expect that global existence result holds only for the dimension d < 4 
when global minimizer for the entropy functional is attained as proven in 
|21j . [25] or [3]. In fact analogous nonexistence results also hold for a problem 
related to fll.ip - fll.4l) but with constant diffusion parameter D (cf. [3]). Thus 
one can conjecture that the gap d G (4,2(1 + v^)) is left for the existence 
of the critical points (possibly unstable) of another type than the extremal 
ones. 

Next, we shall use the aforementioned existence theorems for a given 
temperature 6{t) at time t to prove the existence theorem 14.81 in the mi- 
crocanonical (nonisothermal) setting, i.e. with the given energy and the 
temperature to be determined so that the energy relation fll.l2p is satisfied. 
Steady states for the model were considered, among others, in [22] . Thus, we 
will show that in low dimensions 2 < d < 4 for small mass and domination of 
the thermal energy a gravo-thermal catastrophe (white dwarfs in a physical 
interpretation) does not occur for this system, i.e., neither blow-up for the 
density nor the vanishing of the temperature takes place. 

Finally, in Appendix, we gather the properties of some special functions 
appearing in the Fermi-Dirac model. 

First, notice that due to the self-similar structure of the pressure (11.51) for 
the specific canonical diffusion coefficient D = P' the system (ll.ip - (ll.4l) can 
be transformed to the following one (cf. also the Appendix for Fermi-Dirac 
case and the papers [12], [3], [1] where such D was used). Thus we arrive at 
the system 

nt = V ■ {eP'^ Vn + nP' Vip) in x (0, oo) , (1.7) 

Aif = n in Qx (0, oo) , (1.8) 
{9P'^ Vn + nP' Vv^) ■ u = <^ = on dQx {0,oo), (1.9) 
n(0)=no>0 in C M^ (1.10) 

where we suppose that 2 < d < 4, the temperature is a fixed continuous 
function 9 : [0, oo) [a, b], with some positive numbers a and b with values 
to be determined later. Morover, we look for the solutions of fll.7l) - fll.lOI) 
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satisfying the energy relation given by 

E = ^ [ e'^'^+^P{n9-^/^)dx + ]- [ n^dx= const. (1.11) 

The steady state problem with the prescribed energy for the linear diffu- 
sion was considered, among others, in [2] . The main results of this paper (to 
be specified in the next sections) can be stated as follows 

Theorem 1.1 If P is the Fermi-Dirac pressure and mass M > is suffi- 
ciently small the problem ( [i. 7| )-( I7T7^) admits at least one global weak solution 
for d < 3 and a local one if d = A for a given continuous function 6{t). More- 
over, there exists a local weak solution to the problem ( [i. ?| )-( I7T7^) with the 
energy given by U.ll\) for d < 4. 

Now we shall sketch the method of proving the above theorem. 

First, we regularize the problem to obtain a parabolic system and to apply 
general Amann theory. Next sign-sensitive a priori bounds together with a 
bootstrap argument are used to prove global or local existence depending on 
the dimension d. Then we go with the parameter to infinity and obtain the 
corresponding existence result in weak sense for the original elliptic-parabolic 
system. 

We introduce a new temperature, call it -(9, defined implicitly by the afore- 
mentioned energy relation, i.e. 

E = - [ i^'^/^P{n^-'^/'^)dx + - [ n(pdx= const. (1.12) 

Note the implicit dependence of the 'new' temperature i) on the old one 6 
via n, if (solving (I1.7l) - fll.l0l) for given 9) in the above formula. In sections 
2 and 3, for given 6', we solve fll.7p - fll.10p to get n, yj. Then, in section 4 for 
given value of the energy E, we use the implicit formula fll.l2p for d and ask 
whether the operator T : 6 ^ defined by f ll.l2p has a fixed point. The 
problem of a priori bounds for the temperature 9, determined by fll.llj) . was 
addressed in [6] . In the last section the properties of the special Fermi-Dirac 
pressure function have been gathered. 

Notation. By C we will denote inessential constants, which may vary from 
one line to another. By | ■ |p, for p > 1, we shall denote the standard Lp{Q) 
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norm. By smoothness we shall always mean regularity and it will apply 
only to the function P and is explicitly stated at the beginning of the next 
section. Finally, both || ■ \\hi and || ■ ||* will denote the norm in the Sobolev 
space H^{Q) with and L}{VL) term correspondingly. 



2 The existence result for the perturbation 

In this section we follow the lines of the proof of the existence addressed in 
[3j, where the authors considered a specific Fermi-Dirac density P = Ppo 
defined by (I5.87p . dimension d = 3 and a constant temperature 6, whereas 
here we will just exploit smoothness of the pressure P E {[~^^ oo); [0, oo)) , 
the crucial estimates: a < 9(t) < b with some a > 0, 6 > 0, and for -2 > 0, 

max{po,piz^^'^} < P\z) < ^2(1 + z^''^), (2.13) 
zP"{z) <p^{l + z^/'^). (2.14) 

These assumptions imply that, changing p2 if necessary, 

msx{poZ,piZ^^'^''^} < P{z) < p2{l + 2^+^/^^), (2.15) 

zP'{z) < CP{z), (2.16) 
zP"{z) < CP\z). (2.17) 



In order to study the well-posedness of fll.7p - fll.10p . for A; > 1 and 
P' = P'(n^~'^/^), we consider the following regularized initial-boundary value 
problem 

nt = V ■ {9P''^ Vn + nP' V^) in x (0, 00) , (2.18) 

ift — k A(p = —k n in f2 x (0, 00) , (2.19) 

{eP'^ Vn + nP' Vv?) ■ = = on x (0, 00) , (2.20) 

(n(0),(^(0)) = (no,<^o) in 1) C M^ (2.21) 

with 

no^C^iSl) such that no > , M = |no|i, and (/?o = . (2.22) 
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For this parabolic system we first use the theory developed by Amann [T] 
to prove the local well-posedness of fl2.18l) - fl2.2ip and then the global one. 
The proposition formulated below and its proof has been adapted from [3] 
to cover the case of variable temperature 0{t), slightly more general pressure 
P than the Fermi-Dirac one, and any dimension 2 < d < 4. 

Proposition 2.1 If d > 2 and the function P is smooth and satisfies P' > 
Po > and the temperature is continuous and satisfies 6(t) > a then the 
initial-boundary value problem l^2.18\} - [2.21\) has a unique maximal classical 
solution 

{n,^)eC{nx [0,r^„,);M2)nC2'i(nx (0,T^„,);M2) 

for some Tmax ^ (0, oo]. In addition, n{t) > for t G [Q,Tmax) ■ 

Furthermore, Tmax = oo if there are e > and a locally bounded function 
uj : [0, oo) [0, oo) such that, for every T > 0, the estimate holds 

\\n{t)\\c. + llv^Wllc^ < ^{T) for t G [0,T„,,) H [0,T] . (2.23) 

It is the case for d < ?> if P additionally satisfies Ii2.13\) - ^2.14\ l '^^^ i2.21\) 
with \R'{z)\z^^'^~^^'^ < B and 6{t) < b for some constants b > 0, B > 0. For 
d = A we claim only the local existence result since the boostrap argument 
does not yield \2.23\) . 

Proof. We set Dq = {—6, oo) xM., u = (n, ip) with uq = {uq, (^q), and by the 
assumptions, define a G C^{Do] M2{^)) and / G C'^{Do;R'^) by 

e{P'{ne-'^/^)Y nP'{ne-'^/^)\ f( \=( 



ain) = ^ " ^' -I " -I ' ) . f(u) 

Next, for f G -Dq, we introduce the operators 

d d 

A{v)u = -^^^^diiaijiv) dju) , 
i=i j=i 

d d 

B{v)u = b '^'^h'i {aij{v) dju) + {I2 - b) u , 
i=i j=i 
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where aij{v) = a{v) Sij, I < i, j < d, and 

^^=(J ?) ' ^^(o o) • 

Then, an abstract formulation of fl2.18p - fl2.2ip reads 

Ut + A{u)u = f{u), (2.24) 
B{u)u = 0, (2.25) 
m(0) = Uq. (2.26) 

Thanks to the strict positivity of P'{z) > Pq and the lower bound for the 
temperature 6{t) > a, the eigenvalues of the matrix a{v) are positive for each 
V & Dq, and the boundary- value operator {A, B) is of separated divergence 
form and is normally elliptic in the sense of pj, Section 4] Therefore we 
may apply [H Theorem 14.4 and Theorem 14.6] to conclude that, for some 
Tmax £ (0, oo], fl2.24p - fl2.26p has a unique maximal classical solution 

u = {n,^)eCiQx [0,T^,,);Do)nC2'i(nx (0,^„,,);Do)• 
Also, since no > and the first component of f{u) is equal to zero, the 
comparison principle (see, e.g., [H Theorem 15.1] or [191 Corollary 1.2.1]) 
implies that n{t) > for t G [0,Tmax)- Furthermore, since / does not depend 
on Vm and n > 0, Theorem 15.3 in [T] ensures that T^ax = cxd if there are 
e > and a locally bounded function u : [0, oo) [0, oo) such that (12.231) 
holds true for every T > 0. The fact that the assumptions imposed on P 
guarantee fl2.23p requires some preparatory lemmas and is postponed until 
the end of this section. □ 

We proceed to present a series of lemmas which will guarantee that (12.230 
is satisfied for 2 < < 3 and thus ascertain the global solvability of the 
perturbed problem. We recall after [3] that the neg-entropy functional W 

W = j (^H-{^ + l^ P^"/2^ dx 

plays the role of a Lyapunov functional for the original and regularized prob- 
lem. The function H{z) depending on z = n9~'^^'^ is a primitive of P'{z)/z. 
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However, in our case this functional is not useful for a priori estimate of the 
density n (contrary to isothermal case [3]) as due to (12.271) it is of too low 
order in n. Indeed, the order is 1 — 2/d for the Fermi-Dirac case to be exact 
(cf. [H Lemma 3.6]), and as such does not provide any reasonable a priori 
estimates for the density n. On the other hand, it can be used to get a priori 
bounds for the fixed points of the temperature operator T as was done in [S] 
and is presented in section 4. In the isothermal {6 = const) crucial 
L^+^/'^ bound was obtained from the fact that the entropy (other than W) 
was coercive in this space. As one can see it is not the case for W. For the 
details of the nontrivial derivation of the entropy W one can see [5] and for 
its application to get a priori bounds for the temperature - [6] and [21] . 

Now, we are going to formulate analogous results to the ones presented in 
[21j where a priori bounds for the limit parabolic-elliptic system, as — >■ cxo, 
were obtained. 



Lemma 2.2 Assume that, for d > 2, 

P(z) =piz'+^/'' + R{z), (2.27) 

where the lower order term satisfies \R'{z)\z^^'^^^^'^ < B. Then, for any fixed 
T > and any t G [0,T] fl ^,Tmax), the following growth condition holds 



d / d 
dt V2 



— I - / pin^+^Z"^ + / n<^ 



(2.28) 



Proof. Let t E [0,T] n [0,T^ax) and recall that both P and P' are the 
functions of n9~'^^'^. Now, we multiply (I2.18P by |?7,^/'^ and integrate over 
to obtain 



rf2 d 



Jq Jn Jn 



2{d + 2)dt 

Similarly, multiplying (I2.18P by Aip, we get 



A / ntipdx = -A / P' n \V(p\'^ dx - AO / P'^ Vn-V^dx. 
Jn Jn Jn 
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Summing up the above equalities and using the Holder inequality 
Jn 

taking A{d\ 2)pi = d and 



we arrive at 



-7- ( - /" [ nipdx + - [ iVipl'^ dx] + \ ! ipf dx 

dt \2 Jn 2 Jn J k J^ 

< J_ /■ |V<^|2 (AP'^^V2^-l/rf+l/2 _ ^-l/2^1/d+l/2)2 

4^ Jn 

This yields the claim, with C = -^j^y^, from the assumption on R' apphed 
to the differentiated pressure P'{z) = {pi{d + 2)/d)z^/'^ + R'{z). □ 

Remark. Note that the above theorem holds both in the polytropic case 
with R{z) = and, less obviously, in the Fermi-Dirac case as explained 
below. Indeed, by the properties of Fermi functions (cf. Lemmma 15.61 from 
the Appendix or for more properties see [H Sec. 5]) we get |-R'(-2)| < Bz'"^^^ at 
z = oo and |-R'(-2)| < B at z = 0, implying the required estimate if 2 < d < 6. 

Next lemmas will allow us to estimate the right hand side of (12.281) . 
Lemma 2.3 For any 2 < d < A we have 

[ mpdx <CM^I^-^I'' ( f n^+^'Ux+ [ \V^\^dx\ (2.29) 
Jn \Jn Jn J 

Proof. The proof of fl2.29p involves standard Holder and Sobolev-Gagliardo- 
Nirenberg inequalities as follows 

d+2 



n 



l/2-l/d^l/2+l/d^^^ 



n 



< M^/'^-^/'^ ( / dx 



n 



2d 
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due to the inequality < and the fact that i^o(fi) can be imbedded 
in L^^{Q). The proof of the case d = 2 is straightforward by the Poincare 
inequahty. □ 

In low dimensions a similar argument leads to another estimates (cf. 



Lemma 2.4 For d = 2 we have 



/ rup dx 




Jn 


Jn 



while for d = 3 the estimate 



I 

Jn 



rup dx 



<CM'/' + ^ I n^'/^dx. 



(2.30) 



(2.31) 



holds. 



Now, we are ready to deduce the following lemma on a priori estimates. 

Lemma 2.5 Assume that 2 < d < 4, condition l[2.21\) holds for a smooth 
function P satisfying \R'{z)\ < Bz^^^^^^"^ and the temperature is bounded 
from above 6{t) < b. Then for any t G [0, T] fl [0, T^ax) (^i^d sufficiently small 
data, i.e. mass M if 2 < d < A or the Poincare constant for d = 2, we have 



Pin^+^/Ux+ [ n<pdx + - [ \V^\Ux+l [ [ ^^^dxds<C. (2.32) 
in 2 k Jq Jq 



Moreover, each of the terms appearing on the left hand side of the above 
inequality is bounded and the constant C may depend on the initial data. If 
d = 3 the assumption on the smallness of M can be relaxed due to Lemma 



Proof. Starting with the direct consequence of Lemma 12.31 true for suffi- 
ciently small mass M and large C (if = 2, instead of making mass M small, 
we have to assume that the constant from the Poincare inequality is smaller 
than 2), 

f |V(^prfa; < C f- / pin^^^'Ux+ [ n^dx + - [ dx] 
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we plug this into (I2.28P and integrate with respect to time to arrive, with 
possibly a larger C dependent on the initial value of the right hand side of 
the above inequality, at 

- / pin^+^''^dx+ [ n^dx + - [ \Vip\^dx<Cexp( [ Ce{sy/^ds^ 

which together with an upper bound on 9 ends the proof. Note that (12.291) 
from Lemma 12.31 shows that the negative term mp dx is dominated by the 
positive ones and thus the last claim of the lemma is ascertained. □ 

The integral version of the estimate (I2.28P from Lemma [2.21 follows by an 
argument similar to the one in the proof of Lemma 12.51 and reads 

V{t) < \/(0)exp {c e{sY'^d^ , (2.33) 

where V{t) = f /^pin^+^Z'^ + ^^nip + \^^ | V(/?p. Note that for = 3 we can 
add CM"^/^ in the definition of V and thus relax the assumption on smallness 
of mass M. 

Now, we state similarly as in [3], where only three-dimensional case was 
treated, a lemma on the improved regularity of Vv?- 

Lemma 2.6 Let q,a E (l,oo), d > 2 and T > 0. There is a constant C 
depending on q, a, d and T such that, for t G [0, T^ax-) H [0, T], 

f iVifis)^ ds < C\n\'' ^ . (2.34) 

Jo L'J(0,i;Ld+"{n)) 

Proof. We infer from ^17^ Corollaire 1.1], as in j3j where for a = ^2 '!^% ^^"^ 
d = 3 the authors used the bound with da/{d + a) = 1 + 2/d norm of n, that 

T IV^tl -ikL. + |A09| _da_ < C \n\ _da_ 

k Ll(0,t;Ld+a(^Q)) 'L1{0,t;L'i+a(Q)) L9(0,t;Ld+a (Q)) 

Now, if t G [0, Tmax) n [0, T], we get from the above inequality 



Jo W 'd+c 



ds < C |n|^ 

(C) Li(0,t;Ld+a[n)) 



To conclude we use the imbedding of W^'^{n) in W^''^{n). □ 
Furthermore, an L^-estimate is available for n. 
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Lemma 2.7 Let T > 0, 2 < d < 4, P be smooth and satisfy \2.13\) and 
\2.21\) . There are constants c, C > depending on sufficiently small mass 
M, hounds on 6{t) G [a,b] and P and the initial data such that, for t G 
[0,T^ax)n [0,T], 

Ht)\l + c I \Vn{s)\lds<C (2.35) 
Jo 

In fact, a constant C is a function of the integral \n(s)\i^2/d ds, locally 
bounded in t, which can be estimated by constant due to 1^2. 3 2\) . 

Proof. Note that the estimate of the first term in fl2.35p follows from Lemma 
EHif = 2. Let t G [0, Tmax) n [0, T] and multiply (KTE^ by 2n, and integrate 
over Q to obtain 

^\n{t)\l + 2e / |Vn|2 = -2 [ nP' Vn ■ Vip dx . 
dt Jfi Jfi 

Next, we have 

2 f nP' Vn-Vip dx <e [ P'^ \Vn\'^ dx + ^ / |Vy?|^ rfx 

by the Young inequality, whence 

^\n{t)\l + e [ P'^ \Vn\^ dx<]- [n^lVifil^dx. (2.36) 
dt Jq Jq 

For = 2, by Lemma [2.61 and Lemma [2.51 we deduce that 

|V<^U < C (2.37) 

and that the integrated with respect to the time variable the right hand side 
of fl2.36p is bounded. Thus the estimate fl2.35p is proved in this case. If d > 3 
a longer argument is required. Namely, it follows from the Holder and Young 
inequalities that for any e > 0,a > 2 and some C = Ce 

_„2 a(d+6) 2a(d+6) 

n 6 ' ^ ' 

Then, interpolating with positive {3 = ^^l^^^^y^? we get 

a(d+6) 

< M'\n'^"%„^,.,) (2.38) 
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and P{z) > piz^'^'^/'^, or precisely 6'^/'^~^^ P{n6 '^/^) > implies that 

r p 2c»(d+6) 

j^n' iV^r dx < -M^\e'/'^'P\l,/^,.2) + ClVv^l (2.39) 

1 2a(d+6) 

< -M^||^'='/2+ip||2 ^ ^ , (2.40) 

where the last inequality follows from the continuous imbedding of H^{Q) in 
^2d/{d-2)^^'j constant and the norm 

|2 



||2||^ = y |V;zrc/x+(^y \z\dxj . (2.41) 

Consequently, by fl2.36p . it follows that, for c = 1 — |iVf^, 

rl C C. f . 2Q(d+6) \ 

^|n(t)|2 + c y ^P'2|y^|2 < ^ n^'^/2+ip|2 + iVy^U^"-^'* J . (2.42) 

Now, integration with respect to time, the assumption on growth of P and 
Lemma ED with q = a = 1 + d/2, /? = 0, for 2 < c/ < (3 + Vl7)/2 

(including d = 3) yields the estimate by the time integral of |?T.|i+2/d. Finally, 
due to Lemma 12.51 providing a bound for |n|i+2/d, the estimate (I2.35p is 
proven for any mass M. 

Now, allowing higher dimensions 2 < d < 2(1 + \fi) (including c? = 4) 
we use q = a = ^ > 2, /? = -i^^, da/{d + a) = 1 + 2/d, to 

get the estimate by the time integral of |n|i_|_2/(i but this time for small mass 
only. To get a bound for |ra|i+2/d, assumptions have to be more restrictive, 
e.g. Lemma [2.51 requires 2 < < 4. 

In fact, we have obtained the estimate 

\n{t)\l + c [ eP'^\Vn{s)\l ds < C , (2.43) 
Jo 

which due to the estimate P' > po implies f l2.35p . □ 

Proof of the global existence part of Proposition 12.11 We are now 

ready to prove f l2.23p and thus obtain the global existence. Let T > and 
t e [0,Tjnax) n [0,T]. We claim that there is C > depending on uq and T 
and bounds on 9 such that 

|nP'|^2{d+4)/(d+2)(f^x{o,i)) + I V(riP')U2(nx(o,t)) < C . (2.44) 
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Indeed, we infer from assumption fl^TT^ - d^TTD that zP'{z) < C (1 + ^1+2/^) 
and zP"{z) < C (1 + 2;^/'^) for 2; > 0. Consequently, 

sup \n{s)P'{s)\^+ [ \\n{s)P'{s)\\l^ ds <C (2.45) 

sG[0,t] Jo 

holds. Next, we use the continuity of the imbedding of H^(^l) in L?''^/'^'^~'^\VL) 
and an interpolation argument to deduce f l2.44p . 

We now employ a bootstrap argument to show that (12.231) holds true. It 
follows from (12.151) and the Sobolev imbedding that 

which, together with (12.351) . leads to 

' u Jo 

Therefore, 

\^\L^+S/d{Qx{0,t)) < C , 

and we infer from [191 Theorem IV. 9.1 and Lemma II. 3. 3] that 

i2{d+2)(d+4)/{d2-8)(Qx(o,t)) 

This estimate and (12.441) ensure that 

|V(nP') ■ Vv?|^(d+4){d+2)/(d2+3d)(f^^(o^^)) + \nP' Av9|^(d+4)/(d+2)(f^x(o,t)) < • 
Since 

nt - (P'^Vn) = V(nP') ■ Wif + nP' Aip , (2.46) 
we use once more |19[ Theorem IV. 9.1] to obtain that 

(d+2)(d+4) 

which, in turn, implies that n G L d^+d-e (f2 x (0,t)). With thus improved 
n we would like to bootstrap once again. The right hand side of (I2.46P is in 
the space L'^{Q x (0,t)), with 

q = d{d + 2){d + 4) min {l/(c/3 + Ad^ - 12), l/(2(rf + l)(c/ + 3)(c/ - 2))} . 

14 



Therefore, for = 3 we finally get the right hand side of fl2.46p in L'^{Qx (0, t)) 
with q = 35/12 larger than critical 1 + 3/2 allowing to conclude with 

\\n\yim)<C, (2.47) 

for some e > hj [TS| Lemma II. 3. 3]. However, for d = 4 one should note 
that we have obatined from the boostrap the integrability of the right hand 
side of f l2.46p of the order q = 48/35 which is less than we had before, i.e. 
q = 8/5. Thus for = 4 we cannot conclude with the estimate fl2.47p . □ 



3 The local existence result for the original 
elliptic— parabolic problem 

In this section we shall subtract a convergent subsequence of solutions to 
fl2.18p -( |2.2ip obtained in the previous section which will guarantee the fol- 
lowing existence result for the limiting problem fll.ip - fll.4p as k —>■ oo. 

Theorem 3.1 Assume that M is small enough if necessary, and P is smooth 
and satisfies ^MIW-(KW ^'^^ Moreover, let 2 < d < A, uq e L'^ip) 



and, for given constants < a < h, 9 G C(0, T; [a, 6]). Then there exist 
a weak local-in-time solution n G C (0, T; L^(fi)) , ip G (0, T; if^(i7)) , 
Qd/2+ip ^ (^Q^T;H\n)) of the system (n[j- [LTm . i.e. 

[{n-no)xdx + [ f Vx ■ {OP''^Vn + nP'Vip) dxds = 0, (3.48) 
Jn Jo Jn 

Aif = n, (p = on dn, (3.49) 

for each test function x G H^^'^°'/('^~^)(r2). Additionally, 

\nit)U+2/d + ymm<C, (3.50) 

\nit)\2+ / \\e''^'+\s)Pis)f^,ds<C, (3.51) 



for any t G [0,T) where C depends on uq, VL, d, a and h. If d < 3 then the 
global result can be claimed. 
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Proof. We follow the lines of the proof from [3] , where d = 3 and a constant 
temperature 9 were assumed. We consider uq G L^(fi) such uq > a.e. 
in Q and put M = \no\i (sufficiently small if necessary). Let (no,fc)fc>i be 
a sequence of nonnegative functions in C°°{Cl) approximating ^o, i.e., 

|no,A:|i = M and lim \no^k — ""-ob = . (3.52) 

For /c > 1, we denote by {rik, (pk) the unique classical solution to fl2.18l) - fl2.21l) 
with initial datum (no,fc,0) given by Theorem 12.11 and let Pk = P{nk9^'^^'^). 
Owing to (ES]), fl2:i3D and fl332D there is C > such that 

\nk\2 + yk\\m + l r \{Vk)s{s)\lds+ f \\d'"^+\s)Pk{s)\\l,ds<C. 
I"^ Jo Jo 

(3.53) 

Observe that the Holder inequality, fl3.53p and assumptions (12.131) . (12.161) 
imply 

\OPL^Vnk\^ < C (\n''' VPfcl^ + e''^^+^\V Pkl^A 

d+2 y I /V ^_|_2 d+2 / 

<c(\nu\2\VPk\2 + \VPk\2A , 
\nkPiy^k\d/2 < C\Pk9''/^ Vipk\d/2 < C|Pfc|^|V<^fc|^, 

whence, by (I3.53P thanks to the imbedding of H^(n) in Ld^(Q), 

f (|^(3)Pf (s)Vnfc(s)|2,/(,^2) + \nk{s)P',{s) VMs)\l/2) ds<C. (3.54) 
Jo 

For d = 2 due to (12.371) we get the L°° bound for Vipk whence In^P^Vv^fcli < 
c|Pfc|i < C. We then deduce from the above inequality and equation (12.181) 
that 

l(^fc)i|L2(0,T;tyi.2d/(d-2)(Q)/) < C . (3.55) 

Consequently, owing to (I3.55p . (I2.35P and (I3.53P the sequence (n^) is bounded 
in L'^{0,T;H\n)) and in H\0,T;W^'^'^/'^'^-^\n)'). Owing to the compact- 
ness of the imbedding of H^{Q) in L^(fi) and to the continuity of the imbed- 
ding of L2(^]) in Vri'2««/(^-2)(^])', we infer from [20l Corollary 4] that (n^) is 
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relatively compact in x (0,T)). Therefore, there are n G x (0,T)) 
and a subsequence of {rik) such that rik n a.e. and 

rifc— in L\n x {0,T)) f] C{[0,T];W^'^''/'-''-^\ny) . (3.56) 
Let V2 e L°°(0,T; be the solution to 

A^ = n in ^]x(0,T), ^ = on 9fi x (0, T) . (3.57) 
It follows from (I2.19P and (I3.57P that ipk — ^ solves the Poisson equation 

-A((pfc - y}) = n - rifc - (V9fc)t 

with the homogeneous Dirichlet boundary conditions, and the right-hand 
side of the above equation converges to zero in L^(f2 x (0,T)) by (13.531) and 
(I3.56p . Therefore, 

ypfc — in L'^{0,T;H^{n)) . (3.58) 

Combining (13.530 with the convergence results (13.560 and (13.581) finally allow 
us to conclude that P^^Vn^ and Pj^UkVipk converge weakly to P'^V n and 
nP'V<^ in L2'^/('^+2)(^ X (0,T)) and L'^/^{n x (0,T)), respectively. It is now 
straightforward to pass to the limit as A; ^ oo and conclude that {n, (p) is a 
weak solution to (ll.7p - (ll.9l) as stated in Theorem 13. 1[ 

We may also pass to the limit in (13.500 and use classical lower semicon- 
tinuity argument to deduce that (13.501) holds true. 

Next, by (12.131) and (12.151) it follows from the conservation of mass, (12.420 
and the Poincare inequality that 

rJ , f 2a(d+6) \ 

-|n,(t)|^ + 7 (|n,|^ + ||^'^/2+ip,||^,) <C M +|V^,U--'' j (3.59) 
for some positive constant 7. Integrating with respect to time, we get 

/"* / 2a(d+6) \ 

W{t)\l<\n^A\^~^''rC \ l + |V</^fc(s)U«"-^ e^(^-*)rfs (3.60) 
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for if: > 0. Next, from the Fubini theorem and the double integration of fl3.59l) 
we obtain, for t > 1, 

\\9'/'+\s)Pkis)\\l, ds < j'^ ^ \\d'"^+\s)Pk{s)\\l, dsdr 

1 + Mr)\ldT + J^ \VMs)\a''-' ds\. (3.61) 

Now, I Vv^fcU is bounded in L'^(0, t) for any q G (1, oo) by Lemma [231 and we 
infer from (I3.58P and the continuous imbedding of H'^{Q) in W^'°'{fl) that 
|V(/?fc — V(/5|a converges to zero in L^(0,t). Consequently, by interpolation, 

2a(d+6) 

I V(/9fc — Vv5 1 o converges to zero in L e^-d (0,t). Then one can pass to the limit 
as — s> oo in ( 13.601) and ( 13.611) with the help of ( I3.56P and weak convergence 
arguments for the left-hand sides and conclude that 

/"* / 2a(d+6) \ 

\n{t)\l < \no\l e-^* + C / 1 + \V ipis)\a"'-'' e^^^"*) ds 



for t > 0, while for t > 1 

ft+l / ft rt+2 2c(d+6) \ 

ms)^P{s)\\l^ds<C (1 + Hs)\lds + \V^{s)U'-' dsj 



d{d+2) 
cP-d-2 



Since 93 is a solution to equation ( 11.81) . by ( 13.501) we get taking a 

iVv^U < C ||v9||vi/2,i+2/d < C \n\i+2/d < C. 

Inserting this estimate in the previous two inequalities yields the boundedness 
of |n(t)|2 with respect to time, and then (I3.5ip . For ci = 2 we have 1 + 2/d = 2 
whence the estimate for |n| 1+2/d is sufficient. □ 



4 Fixed point for the temperature operator 

r 

First we recall a lemma on relations between n and imposed by (I1.12p . 
This should be understood as necessary condition for the density obtained 
from ( ll.7p -( fLT0l) and not as a sufficient condition for admissibility of the 
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given energy E. The lemma on a priori bounds is related to the one from 
[3] in the Fermi-Dirac case and to the ones from [5J and |6] in more general 
case. Recall from [6], Lemma 3.1] or [21] the following version of these energy 
estimates. 

Lemma 4.1 Let v = 4/(ci(4 — d)) for any 2 < d < A. Provided that P{s) > 
Pis^^'^^'^ for any ^pi > e > and all s > 0, the following estimate holds 



E + CM^^^ > max < e n 



(4.62) 



Moreover, for each < e < d/2, the temperature and the density n should 
satisfy 



E>e[ ¥^'^+^P{n^-^l'^) dx+ f dx 
Jn Jn 



CM 



(4.63) 



Now we shall prove some a priori estimate for norm of the solution 

to BVP fll.7l) - fll.lOI) . We derive them directly from these equations since at 
this moment we cannot directly use the energy a priori bounds presented 
above. Note that these a priori estimates for limit functions are better than 
those for the perturbed parabolic system presented in previous sections (cf. 
lemmas: [221 [231 [23]) . 

The next lemma can be found in [21| (cf. Lemma 2.1 and 2.2 therein). 

Lemma 4.2 For any 2 < d < 2(1 + \/2) and (p related to n by ( fi.gj) we have 
the estimate 

(4.64) 



(4.65) 



/ n(p dx 


< CM^'^/'^ [ 


Jn 


Jn 



Let 2 < d < 4 and assume that 

P{z) =piz'+^/'' + R{z), 



where the lower order term satisfies \R'{z)\z^^'^ ^/'^ < B. Define the 'asymp- 
totic energy', i.e. E°'{t) = lim^^o+ -E'('^); 



E\t) 



pi n^+^''^ dx 



\V^\^dx. 



(4.66) 
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Then for any fixed T > and any t G [0, T] the following growth condition 
for E"" is available 

^E'^it) < Ceity/^ [ \Vip\^ dx. (4.67) 
dt 

Remark. Note that the above theorem holds both in the polytropic case 
with R{z) = and in the Fermi-Dirac case, since by the properties of Fermi 
function (cf. Lemma [5.61 the Appendix and [H Sec. 5]) |-R'(2;)| < Bz^'^/'^ at 
z = oo and |-R'(-2)| < -B at z = 0. 

Remark. It should be noted that for the polytropic case the theorem implies 
the dissipation of the energy, since in this case E"' = E. 

Applying the estimate (14.641) to mp dx = — I Vy^p dx and integrating 
(I4.67P from Lemma l42l allows us to derive the following corollary (for details 
and the proof see 



Corollary 4.3 Under the assumptions of Lemma \4 ■ S\ E°' grows like 

E%t) < E%0) exp {C{t)) (4.68) 

where the function C is defined by C{t) = cM'^^'^/'^ j^{0{s)Y/'^ds. Moreover, 
E"' is positive if CM^~'^^'^ < dpi while for d = 2 we assume smallness of the 
Poincare constant i.e. C < 2pi. 

After integrating inequality (14.681) from Corollary 14.31 and using Lemma 
14.11 we obtain L^"*"^/"^ estimate for the density n. 



Corollary 4.4 Under the assumptions of Lemma 4-^^ we have for any 2 < 
d < 4, u = 4/((i(4 - d)) and any M > 

[ n^+^^Ux < CMi+^ + |E"(0)|exp (C(t)), 
Jn 

while for small M > and any 2 < d < 4 the constant C may depend on M 

n^+^/rf dx < CE"(0) exp (C(t)) 

n 

where CM^"^/*^ < dpi . 
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Corollaries 14.31 and 14.41 allow us to to define the new temperature This 
was the subject of the considerations in [21] under physically acceptable 
property of the pressure || > expressed as 

P{z)z-^-^^'^ \P2>0 (4.69) 

that guarantees, in particular, the uniqueness of the temperature emerging 
from the energy formula (11. lip . 

The next theorem claims that that temperature is well defined for some 
values of the energy (for the proof see Theorem 3.2 in [21]), and the remainder 
of the section is devoted to proving its compactness. 



Theorem 4.5 Assume that P is smooth and satisfies ^2.13\) . \2.14\j , \4-69{ ) 



and 1^2.21^ Then the temperature operator T : ^ d is formally well defined 



by U.ll\) for small mass M <C 1 and 2 < d < A all the values of the energy 
E admissible at t = 0. Moreover, for 2 < d < 4, u = A/ {d{A — d)) and some 
positive constants B,C, it has to satisfy 

E < BM^+^I'^ - CM^+''. (4.70) 

Next, we estimate to get the compactness of the operator T. By 
differentiation of the energy relation (11.111) we get 

^'K/J""'^^^^I"*"'^)(/J''f' ■ ''''' 

In the following two lemmas we claim the boundedness of both factors in 
appropriate norms so that be in L'^ with some 7 > 1 which guarantees the 
equicontinuity condition in the classical Arzela-Ascoli thoerem. 

Lemma 4.6 Assume that P is a smooth function such that 

{P{z)z-'-^^''y <0, (4.72) 
{{P{z)z-'-^/''yz'+'/'')' > 0. (4.73) 

Then the first inequality implies Ii2.16\) with a strict inequality and C = 
1 + 2/d, i.e., 

P\z)z < (1 + 2/d)P{z) . (4.74) 
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Moreover, the function p is decreasing, convex with respect to and satisfies 

%dx>C (4.75) 

for some C > depending on M and a lower bound for d > a provided that 
liminf (-P(2)2-i~2/dy > > 0, (4.76) 

z^oo 

liminf {-P{z)z-^-^/'^)' z^''^+^ > C > 0. (4.77) 
Proof. The formula for the first derivative reads 

g = ~{dl2)r" {P'{z)z - (1 + 2/d)P{z)) , (4.78) 
or in an another form 

^ = -(rf/2)?9^/2^2+2/d (p(^)^-i-2M' > , (4.79) 
ov 

where z = n'd~'^^'^. Then, the second derivative can be calculated 

{2/dfd^-'"^-^ = P"{z)z^ - (1 + 2/d)P'{z)z + (1 + 2/d)P{z) , (4.80) 
or expressing it in a more concise way 

(2/rf)2^9i-^/2^ = z-" ((P(^)z-i-2/'^)'^i+2/dy _ 

Thus the convexity of p with respect to follows from the second assump- 
tion (which by the way can be deduced from the first assumption or (14.741) 
under an extra convexity assumption on P). Now, by the asymptotics of P, 
i.e. gZED, dTTZD, liminf J^^dx> CM as and lim inf ^ || ^a; > 
CM^~^/'^a, respectively. Hence, by the convexity of p with respect to 'd, 
(|4T5|) follows. □ 

To get the bound for i)' we are left to estimate the denominator in fl4.7ip . 

Lemma 4.7 The denominator appearing in (^^77) is bounded in some 
with 7 > 1, i.e.. 



T 



dp 1 d 

—ntdx + -— I mpdx 
on d dt 



7 



dt<C . 



22 



Proof. First, recall that 



Using (13.551) we get the bound for ut in L^(0, T; W'^^^'^'^^^'^^'^^fiy) so we have to 
show that both ^ and are bounded in L^{0,T; with some 



^ > 2. Using (12^ and (13.511) with a = and any g > 1 we get 

f\'^^{sW_^ds<C f\n{s)\lds <C'. 
Jo ^-'^ Jo 

whence cp G L«(0, T; with ^ > 2. Moreover, ^ = OP' ^ 

n^l'^ and the function = Q^^'^l'^p" is bounded in view of the regularity 
assumption on P, so the claim is guaranteed by the estimates for P (12.131) 
and flCTj) and Lastly, from (|232|) follows the bound for l^^\dx and 

from fl2:29|) and Lemma [23] for |V(^pdx. □ 
Finally, we recall after [6J and |i24j a 'priori bounds on the fixed points 
of the compact operator T thus guaranteeing the existence result for the 
problem. 

The authors assumed therein, for negative initial values of the entropy 
liminf {H{z) - {d/2 + l)P{z)/z) > W(0)/M , (4.81) 

where H'{z)z = P'{z) under the following conditions consistent with (12.131) 

P{z)/z^+^/'' \ £ > 0, (4.82) 
liminf P(2)/z > 0. (4.83) 

z\0 

If (I2.27P is fulfilled then the highest order terms cancel assuming that the 
limit exists 

, Hiz)z- {d/2 + l)P{z) , / , d ,,,, ,\ 

lim = lim H{z) - -zH\z) (4.84) 

= lim {G{z) - {d/2)zG'{z)) = Gq (4.85) 

z^oo 

where 

H'{z) = hi{2/d)z^/''-^ + G\z) (4.86) 
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with G'{z) = o[^z^l'^ we are left with the analysis of the lower order term 
G\z) = gizf^ + o(z^) with some (3 < 2/d ~ 1. Namely, if (3 < —1 then 
Go = in fl4.84p . e.g. for the Fermi-Dirac model in d = 2. Otherwise, if 
j3 G [—1,2/d — 1) then the important factor is the sign of gi which has to be 
positive, and indeed is, e.g. gi = 1 — 2/d and P = —2/ dm. the Fermi-Dirac 
case for d > 3, to imply Gq = oo and thus to guarantee (14.81 1) . 

Thus we have proved the existence of a fixed point for temperature op- 
erator T and we can formulate the following existence result in the micro- 
canonical case. 



Theorem 4.8 Assume that P is smooth, satisfies ^2.13\) . i2.14\) , i2.21\) , 



\4-69\} , U-82\ ), U-83\ ) and U-84\ ) '"'^^^ Go>0. Then for the negative initial 



values of the entropy we get the global existence result for ( 1.7 )- fTT^) with 
the energy constraint 

5 Appendix on Fermi-Dirac model 

First, it should be noted that for the Fermi-Dirac case we have 

dPpniz) = (/,-/2_i(2;.//i)) , (5.87) 

where fx = T]oGad2'^^'^, G is the gravitational constant, r/o ^ a bound for the 
density in phase space and fa is the Fermi function of order a > —1 defined 
by 



Uz)=J^ Y^^dx. (5.88) 

In [H Lemma 5.1] substitute z = — log(A) and fa{z) = Ia{e~^) to get 
Lemma 5.1 The following asymptotic relations hold as z —>■ oo 

fa{z)-^ = 0{z^''), (5.89) 
for each a > 0, while for each a > — 1 

z-'^ lfaiz)z - ^fa+iiz)] - . (5.90) 
a + 1 3 
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Moreover, we have the recursive relation for the derivatives 



fLiz) = afa-i{z). 



(5.91) 



Next P, Lemma 2.2] can be formulated as follows. 

Lemma 5.2 For a > (3, fa o f^^ is an increasing convex function. 

In conclusion of the above lemma the function Ppo shares the same prop- 
erties. 



Next, to check that the assumptions of the Lemma are verified for the 
Fermi-Dirac, case we will need a version of [H Lemma 5.3]. 

Lemma 5.4 For all (3 < a + 1 the following inequality holds 



Now, we recall after [3, Lemma 2.1] the following properties of the Fermi- 
Dirac pressure. 

Lemma 5.5 The function Pp£, belongs to C^([0, oo)), is nonnegative, in- 
creasing, convex and can he extended to an element (still denoted by Pfd) of 
C^{[—6, oo)) for some 6 > 0. 

Next the following asymptotic result holds (cf. [21]) for 



Lemma 5.3 The function Pfd, defined by i\5.87\ ), is increasing and convex 



function. 




(5.92) 




Lemma 5.6 For the Fermi-Dirac pressure Ppo we have at z = oo 




(5.93) 



and, in consequence. 



where pi = -^{d/ jj)'^^'^. Moreover, at z = 0, we have 



(5.94) 




(5.95) 
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Lemma 5.7 The conditions || > and/or P{z) j z^^'^l'^ \ P2 > are sat- 
isfied for 

Pfd{z) = (/i/rf)/d/2 o f-/l_, {2z/fi) . (5.96) 
Proof. Indeed, putting 2z = ^fd/2-i{x), we get 

^ = ((-/^/2)(/./2-i(a:))^(/^/2_i(x))-^ + (1 + 2/d){^^/d)f,/,{x)) . 

Thus, the condition || > is equivalent to 

- ((/x/2)/,/2_i(x))' + (/i/rf)/^/2_i(x)(/x/2)(l + 2/rf)/,/2(x) > . 



This, however, follows from Lemma 15.21 (take a = d/2 and (3 = d/2 — 1) 
or more explicitly by the property of Fermi functions presented in fl5.92p . 
namely ^ U^) < . □ 



Lastly, we shall trace how flL7l) -( |TTT0l) could be derived, in the Fermi- 
Dirac case, from flLll) - fll.4p (used by the authors in [3]) under the assumption 
(11.51) with a specific diffusion coefficient, used in [3] and [12], 

^(^' = TF^' (5.97) 

where Ia{e~^) = fa{z) and A = ■ Using the recurrence property 

(I5.9ip (cf. also |4j, Section 5] and |i21j, Lemma 1.1]), we get, differentiating 
formula (15.870 . the relation D = P'. Furthermore, OD = Moreover, it 
should be noted that in [3] and [1] the authors used the following notation 
F' = P'^, V = nP'. Note that D should be defined exactly as in (K97^ but 
it might differ throughout these papers up to a constant, inessential therein. 
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